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Abstract 

The spatial momentum dependence of the spectral function for vr and a at finite temperature is 
studied by employing the 0(4) linear sigma model and adopting a resummation technique called 
optimized perturbation theory (OPT). The poles of the propagators are also searched for. 

We analyze the spatial momentum dependence of the imaginary part of the self-energy and find 
its temperature-dependent part to vanish in the large momentum limit. This is because the energy 
of the particles in the heat bath which participate in the process becomes large and therefore the 
Bose distribution function vanishes. 

We then calculate the spectral functions and search for the poles of the propagators. First, we 
discuss the temperature dependence for zero spatial momentum. We reproduce the spectral func¬ 
tions in both TT and a channels in the previous work and find that the pole, which is responsible 
for the threshold enhancement of the spectral functions in both tt and a channels, is different from 
the one for the peak at zero temperature. Secondly, we discuss the spatial-momentum dependence. 
When the temperature is low the spatial momentum dependence is also small. As the temperature 
becomes higher the spatial momentum dependence becomes also large. When the spatial momen¬ 
tum is smaller than the temperature, the spectral functions do not deviate much from those at zero 
spatial momentum. Once the spatial momentum becomes comparable with the temperature, the 
deviation becomes considerable. As the momentum becomes much larger than the temperature, 
thermal effects effectively decrease. 
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I. INTRODUCTION 


One of the current interests in the strong interaction physics is how the properties of 
hadrons change at hnite temperature and/or density. To study this has become increasingly 
important since it is closely related to the data obtained in currently carried and future 
planned high energy heavy ion collision experiments [1]. 

A novel feature at hnite temperature and/or density is the absence of Lorentz symmetry 
due to the existence of the heat bath and/or matter and therefore, the spectral functions of 
hadrons depend on two rotationally invariant variables in addition to temperature and/or 
density, for example four momentum squared and spatial momentum squared. Since the 
hadron spectra observed in the heavy ion collision experiments have various spatial momenta, 
it is important to know theoretically how the spectral functions depend on the spatial 
momentum in order to draw a conclusion about the change of the properties of hadrons 
at hnite temperature and/or density. 

The purpose of the present paper is to study the spatial momentum dependence of the 
spectral functions and the pole of the propagator in the complex four momentum square 
plane for hadrons at hnite temperature. We employ a low energy ehective theory of QCD, 
0(4) linear sigma model and adopt a resummation technique called optimized perturbation 
theory (OPT) to it [2]. As is well known, naive perturbative expansion breaks down at 
hnite temperature and proper resummation of higher order terms is necessary [3-5]. By 
calculating the self energies of tt and a in the OPT to one-loop and analyzing their spectral 
functions, we study how the spectral functions depend on the temperature and the spatial 
momentum. We also search for the poles of the propagators. 

The properties of a at hnite temperature at rest were studied in [2, 6-8]. The authors 
in [2, 6] analyzed the spectral function in 0(4) linear sigma model. By adopting OPT 
Chiku and Hatsuda found threshold enhancement of the spectral function due to the partial 
restoration of chiral symmetry. Patkos et ah applied large N expansion, and obtained 
results similar to Chiku and Hatsuda. Dobado et al. studied the mass and width using 
chiral unitary approach. They searched for the poles of the T-matrix of tttt scattering and 
showed that the width grows with temperature, while the mass decreases slightly. 

This paper is organized as follows. In the second section, we give some basics of 0(4) 
linear sigma model and thermal held theory and briehy review the formalism of OPT. Then 
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the origin of the momentum dependence is discussed in the third section. The fourth section 
is devoted to show the numerical results of the analysis. Finally we summarize this paper 
in the fifth section. The explicit expressions of one-loop self energies at finite temperature 
are given in the appendix. 


II. FORMULATION 


A. Linear Sigma Model 

The Euclidean Lagrangian density of the 0(4) linear sigma model is given by 
Ce = ^[(<90)^ + ~ ^00 + (counter term), 


( 1 ) 


where 4>{x) = (00, tt), and /i0o is an explicit chiral symmetry breaking term. When the 0(4) 
symmetry is spontaneously broken, 0o has an expectation value: 


(0o)/3 — 


( 2 ) 


where <>/? is the thermal average, and 0 = ^. After taking out the expectation value, as 
00 = + C) oiie can rewrite the Lagrangian density Eq.(l) as 


-I-(T^TT^ H—-I —+ (-1" h 'C — h)a + —UE H-;- M 

12 6 6 ^ 6 ^ ^ ^ 2 ^ 4! ^ 

-|-(counter term). 

Here, moo- and movr are the tree-level masses: 

2 2 I ^2 2 2 I ^2 

= /^ + 2^ ’ = h + • 

^ is determined by the condition (a)i 3 = 0 which reads in the one-loop order as 

- 1 )) 


( 3 ) 
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where k is the renormalization point, uj^ = \lk‘^ + and n(uj) 


( 5 ) 
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B. Thermal Field Theory 


In the imaginary-time formalism of the thermal held theory [9], the Matsubara propagator 
is dehned by 

A{iujmP,T) = J y ir, a;)0(O, 0 ) >/3, (6) 

where 

Tr(j){-iT, a;)0(O, 0 ) = 6'(r)0(-ir, a;)0(O, 0 ) -h 6*(-r)0(O, O)0(-ir, x), (7) 

and Un is the Matsubara frequency, Un = 2nnT (n is an integer). The Matsubara propagator 
is dehned for discrete pure imaginary energies only but one can dehne the propagator, 
A{pq,p,T), for any continuous complex energy po, as the proper analytic continuation of 
the Matsubara propagator [3]. Then, the spectral function at hnite temperature, p{p,T), is 
obtained as the discontinuity of the propagator: 

ipiP,T) = A{po + ip,p,T) - A{po-ip,p,T). (77 ^ 0+), (8) 


Introducing the self-energy Il{pQ,p,T) by 

Mpo,P,T) = — -j 

pZ _ ^z 

one can write the spectral function as 


-1 

-n(po,P,T)’ 


. ^_ -21mU{po + ip,p,T) _ 

^ — w? — Ren(po + ip, P, T)y + (Imn(po + ip, P, T)y 


( 9 ) 


( 10 ) 


C. OPT 


It is known that naive perturbation breaks down at T 7^ 0, and resummation of higher 
order terms is necessary [3-5]. We adopt here a resummation technique called optimized 
perturbation theory (OPT) [2]. The idea of the OPT consists of dividing the mass parameter 
of the Lagrangian into two pieces, /i^ = — (m^ — /i^) = — y, and treating the one, 

as the tree-level mass term while the other, y, as perturbation. 

The Lagrangian of the linear sigma model Eq.(l) can be rewritten as follows: 

1 ->1 A 

Ce = -[(<90)^ + - -y0^ -F “ ^0o + (counter term). (11) 
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The thermal effective potential is written as the functional integral: 


V{ 0 ^T) 


ln/T>0exp[-|/(f d^x{C{$+0) + J ■ 0)\ 

/(f d^x 


( 12 ) 


Here, J = ^, Jo d'^x = Jq dr f d^x, and the parameter, <5, is introduced for convenience 
(5 = 1). The OPT is dehned by the expansion in <5 with the assignment, = 0(6^), 

X = 0(6) while the naive loop expansion corresponds to the assignment, = 0(6^), 

X = O(50). 

When spontaneous symmetry breaking takes place, tree level masses are 


2 2 I f-2 


2 2 I ,-2 

mo^ = m + 2^ • 


(13) 


The parameter ^ is determined by the stationary condition: 

dV{(p,T, m^) 


d(f 


= 0 . 




(14) 


Note that the derivative with respect to ^ does not act on This condition is equivalent 
to the condition (a) 13 = 0 in Eq.(5). If any Green’s function is calculated in all orders 
in OPT, they should not depend on the arbitrary mass, m. However they depend on the 
arbitrary mass, m, if one truncates perturbation series at certain order. One can determine 
this arbitrary parameter so that the correction terms are as small as possible. We adopt the 
following condition [2]: 


p"-mL-{n".(p) + nno,r)} 


Jp=(mo7r,0) 


= 0 , 


( 15 ) 


where n°(p) and Il^(p,T) are, respectively, T independent part and T dependent part of 
the self-energy for tt. Explicit forms of a and tt self-energy at one-loop level are given in 
Appendix. 


D. Physical Conditions 

The renormalized parameters /i^, A and h can be hxed by three physical conditions at 
zero temperature. We use the following three conditions [2]. 

(l)The on-shell condition for the pion: A~^(p^ = T = 0) = 0 or 'n-,t(p^ = T = 0) = 0, 
where = mo-,, = 140MeV. 
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(2)Partially conserved axial-vector current(PCAC) relation: firml = is the finite 

wave fnnction renormalization constant which depends on the renormalization scheme, 


Z-rr — 


An,(p2 = m;,r = o) 


is the pion decay constant, = 93MeV. 

(3)The peak energy of the spectral fnnction in the sigma channel is determined from the tttt 
scattering experiment to be = (550MeV)^. 

The parameters given by solving these condition are = —(284MeV)^, A = 73.0, h = 
(123MeV)3, K = 255MeV. 

E. Analytic Continuation 


From now on we regard the propagator as a fnnction of p^, \p\ and T, i.e. A(p^, \p\,T). 
To know the location of the poles of the propagator in the complex plane helps us 
to understand the behaviour of the spectral function. For this purpose it is necessary 
to obtain the analytically continued propagator in the unphysical sheet. Therefore, from 
Eq.(9), the problem is reduced to performing analytic continuation of the self-energy into 
the unphysical sheet. The self-energy has branch cuts on the real axis, p^ < {rriQ^ — mo 7 r)^ 
and p^ > (moo- -l- for < 0 and p^ > 4mg^ for \p\,T). 

Analytic continuation of the function through different cuts define different analytically 
continued functions on different Riemann sheets. Here, we discuss analytic continuation 
of IpI)^) through the cut, 0 < < (mo,r — as an example. Our method is 

the same as the one in ref.[7]. In the physical plane one may define n,r(p^, |p|,T) by the 
following /‘^^^(p^, IpI, moTT,nrocr,T) introduced in Appendix, 


|p|,mo,r,mo<7,T) = 


e 


'- + 2-a^ ln(—1^) - ln(—yi) - c In 


(c -f Ott) (c + Oo-) - 


where 


2 r 


_ dkk / n(u;^k) , ((a^p^ + klpjf-p^u;^/,) 

Jo 167r2|p| \ ^ - k\p\y - plulk) 

I njuak) , ((aaP^ + fc|p|)^-PoOxgfc) \ 

^ o;., {{a.p^ - k\p\y - plool,)) ^ 


aa 2 K 
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(18) 


c = 


5\ 


1 - 


(mp^ + mp^ 

p2 


1 - 


ij^Ocr ^Oti 

p2 


This function is analytic in the complex physical plane except for the real axis. Let 
^unphysiv^i \p\-,^o-K-i^OaiT) be the analytic continuation of |p|, mo,r,into the 

unphysical sheet through the cut 0 < < (mp^ — mpo-)^ and F(p^, |p|,mpT^,mpo-,T) be the 

difference of the two functions on different sheets: 


F(p^ |p|,mp^,mo<„T) = lunphysip"^, \p\,mon,moa,T) |p|,mp^,mp^,T). (19) 

On the cut 0 < < (mpo- — mp^r)^, /*^^i(p, mp^, mp^, T) approaching from the above and 

/<2\ 

^unphysiP^^On,moa,T) approachiug from the below are continuous: 

+ iri, \p\,mo^,mo^,T) = I^nphysip"^ -'^V, \p\,mon,moa,T), (20) 


and therefore, 


F(p^ IpI , mp^, moa, T) = (p^ + ip, \p \, mp^, moa, T) - (p^ - ip, \p\, mp^, moa, T) 

= DiscI^‘^\p^,\p\,mon,moa,T) 

= 2ilmI^‘^\p‘^,\p\,mo^,mo„,T) 


In 


(1 




)(1 




87r(3|p| (l-e-l>“^)(l-e-'3“J)' 


( 21 ) 


where 


u}^ =-a-^^Jp^+ \p\‘^ ±c\p\, = a^^Jp^ + \p\^ ±c\p\, (22) 

which can be trivially extended for complex p^. Thus, the complete expression for 
^ulphysip"^^ \p\,mo^,moa,T) reads 


lr>l m m T^ — ^ i 2 n n i n7r)(c + Qq 

^unphysKP 1 iPli nipTr, mp^, J j IOtT^ 6 ^ Ott ^2 > ^ — O ) (c — O 

/•°° dkk / njuJnk) . ((g^p^ + k\p\)^ - plul^,) 

Jo 167r2|p| ^ {{a^p^ - k\p\y - plult,) 

^ njuj^k) ((g^P^ + k\p\Y - pIuI,,) \ 


uJak ((g^p^ - k\p\y - plulk)) 


In 


(1 


_ p /^i^TT 


)(1 


_ p 0^a 


Snplpl (1 - )(1 — ) 


( 23 ) 



III. IMAGINARY PART OF THE SELF-ENERGY 


The momentum dependence of the spectral function is caused by the self-energy n(p, T). 
Since the real part of the self-energy is related to its imaginary part by the dispersion 
relation, the origin of the momentum dependence of the spectral function can be traced 
back to the imaginary part of the self-energy. Therefore, we discuss here the imaginary part 
of the self-energy [10]. We hrst consider the self-energy for the pion held: 

Imn^(p,T) 

= -sign(po)-^ / TTT-^-, -1 [(1 + ^vr)(l + n^) - n^n^]6{po - 

9 J (zTT)'* 4:UJnUJa- '' 

+ [(1 + n^)n^ - (1 + n^)na]S{po + 

+ [(1 + n^)na - (1 + na)n-^]6{po - 

+ [n^na - (1 + n^){^ + na)]S{po + + a;^))}, (24) 

where Utt = n{uTr), = n{ua), uj-k = ^-|- ttiq^ and Ua- = \J{p — ky + The hrst term 
represents the probability for the process, n na, with the statistical weight (l+n^)(l+no-) 
minus the probability for the inverse process, ttct —> tt, with the weight n^Uo-, which is 
nonvanishing when p^ > (moo- + ■ (When we describe processes we use the notation a 

and TT for on-shell particles while a and fr for oh-shell particles associated with sigma and 
pion helds.) The second term is due to the scattering of fr with particles in the heat bath 
and represents the probability for fTTr —a with the weight + no-) minus the probability 
for a ^ TTTT with the weight no-(l + n,r)- It exists only at hnite temperature. If rrioa > movr, 
the process, nn —> a, is possible when p^ < (moo- — "nroTr)^- The third term and the fourth 
term are anti-particle counter parts of the second term and the hrst term, respectively. 

Imaginary part of the self-energy for sigma has a structure similar to that for pion, but 
the relevant processes for the former are a ^ nn and a ^ aa instead of fr —ttu for the 
latter. Hereafter we assume po > 0. Carrying out the integrations in Eq.(24) we obtain 

Imn^(p,T) 
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= < 


A^e^ 

1 

9 

IOtt p , 

A^e^ 

1 

9 

IOtt p , 

A^e^ 

1 

9 

IOtt p 


+ 

A^e^ 

1 


= 


9 


/ _ duj^ [1 + n(a;^) + n{po - o;^)] > {moa + rrio^f 

JUJ^ 

\- 

du^ [n{uj^) - n{po + o;^)] 0 < < (mo^ - mon)"^ 

OO 

_ du!^ [n(u^) - n(po + 

'tt 

^ duj„ [n{u^) - n{po + o;^)]| -(mg^ - < p^ < 0 


p2 < -{ml^ - ml^) 

p2 > (moa + monY 

< p^ < (moa - mon)‘ 
p2 < -(mg, 


OO 

PCSD 

+ / _ da;, [n(a;,) - n(po + a;,)] 

JiOa 


1 j (1 — e ^‘^^)(1 —e 


1 

Stt 167r/3|p| (1 — 

1 , (i-e-/3-^)(l-e-/'-^-) , 2 ' 2 ^ 


-Ill - 

9 167r/3|p| (i-e-d^^)(l-e-d-^-) 

i^2^2 ^ (^X — 

9 167r/3|p| (1 — e“d‘-^^)(l — 


where and are 


OJZ 


|a^Po| ± c|p||, 


a;; 


= lla^pol ± c|p||, 


and, a^, and c are shown in Eq.(18). 

We hrst discuss the behaviour of the imaginary part at large momentum. As \p\ 
the temperature-dependent part of the imaginary part behaves as 

1 


Imn^(p,T) 




9 167r/9|p| 
^ 


(25) 


(26) 


CX), 


(27) 


This can be understood as follows. As |p| —cx), |a^ ± c||p| and a;, —|a, ± c||p| 

therefore the energies of tt and a, which participate in the process, become also large. Since 
the temperature-dependent imaginary part of the self-energy includes a Bose distribution 
factor for either tt or a, which vanishes at the large energy limit, the temperature-dependent 
imaginary part of the self-energy also vanishes when |p| —cx. This implies that thermal 
effects effectively decrease at large momentum. 

It should be noted, however, that the temperature dependence of the self-energy does not 
completely vanish at large momentum: the contribution of the thermal tad pole diagrams 
to the real part of the self-energy and the temperature dependence of the mass parameter 
due to the optimization procedure do not vanish at large momentum. 
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We next ask at which momentum the imaginary part of self-energy becomes maximum. 
The process which generates the imaginary part is nn a. At low temperature, the 
distribution function n(a;) rapidly damps as |p| increases, i.e. the heat bath is occupied by 
pions almost at rest, while tt can interact with only pions whose energy lies between uj~ 
and . This means that the integral in Eq.(25) has a sharp peak at the momentum |p| 
which satishes ci;^ = mo 7 r, i.e. |p| = cp^/mo 7 r- This momentum maximizes the imaginary 
part of the self-energy, which is proportional to the integral times the kinematical factor 
l/|p|. This is because around \p\ = 0 the integral is proportional to |p| and cancels the 
kinematical factor, which makes the imaginary part of the self-energy constant, while around 
Ip I = cp^/moTT the integral has a sharp peak and the dependence of the kinematical factor 
can be neglected. When the temperature is not low, the momentum which maximizes the 
imaginary part cannot be expressed in such a simple form and is determined by balancing 
the integral and the kinematical factor l/|p|. 

A comment on higher loop effects is in order here. Let us consider the process tt —ttct. 
In Eq.(25) the term is nothing but TT* where T is the tree-level T-matrix for the 
process, n na. It is speculated that one can take account of higher loop effects by 
replacing the tree-level T-matrix by the one including higher loop effects. If this is the case, 
the above argument on the momentum dependence of the imaginary part of the self-energy 
holds beyond the one-loop order as far as the T-matrix is not singular at large momentum. 
Of course, new processes such as ttct —ttct exist in higher loop terms, which change the 
structure of the spectral function. This is a future problem. 

IV. NUMERICAL RESULTS 

A. T dependence for p = 0 

In this subsection we discuss how the spectral functions and the poles of the propagators 
depend on the temperature with the spatial momentum kept hxed (p = 0). 

1. TT channel 

Fig.2 shows the temperature dependence in the tt channel. When T = 0 the propagator 
has a pole, (A), at = (145MeV)^ on the physical sheet of the complex plane which shows 
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up as a (5-function peak in the spectral function. Another pole, (B), exists at = (447MeV)^ 
on the unphysical sheet which is too far from the physical region to have a visible effect on 
the spectral function. As the temperature increases, the pole (A) moves down and away from 
the real axis and the peak of the spectral function becomes broader. The pole (B) moves 
left on the real axis and approaches the fTTr —a threshold as the temperature approaches 
T = 165MeV which causes the enhancement of the threshold. When T = 165MeV, the pole 
crosses the nn ^ a threshold at = (moo- — and appears on the physical sheet. 

When T > 165MeV a bound state peak exists. 

2. a channel 

Fig.3 shows the temperature dependence in the a channel. When T = 0 the propagator 
has a pole, (C), at p^ = (556MeV)^ -f f(368MeV)^ on the unphysical sheet of the complex 
p^ plane which shows up as a broad bump around p^ = (550MeV)^ in the spectral function. 
Another pole, (D), exists at p^ = (207MeV)^ on the unphysical sheet of the complex p^ 
plane which seems to provide a little shoulder at the a tttt threshold. As the temperature 
increases the pole (C) moves left and a little bit down and the bump becomes broader and 
shifts towards lower p^. The pole (D) moves right on the real axis and approaches the 
a TTTT threshold as the temperature approaches 145MeV which causes the enhancement 
of the spectral function at the a tttt threshold. When T = 145MeV, the pole crosses the 
d —TTTT threshold at p^ = and appears on the physical sheet. When T > 145MeV a 
bound state peak exists. 

A comment is in order here on the threshold enhancement of the spectral function and 
the behaviour of the poles of the propagator. What is observed above is nothing but the 
well-known phenomena in the quantum mechanical s-wave scattering [11]. 

Now, we would like to compare our results for p = 0 with previous results. In ref. [2] 
Chiku and Hatsuda calculated the spectral functions in the vr and the a channels using the 
OPT. We just reproduced their results for the spectral functions. We also studied the pole 
of the correlation function in the complex p^ plane, which was not done in ref. [2]. In ref. [6], 
Patkos et ah applied large N expansion, and obtained results similar to ref. [2]. Then, in 
ref. [7] they analyzed the relationship between the scalar-isoscalar spectral function and the 
second Riemann sheet pole in the same channel. The behavionr of the pole around the 27r 
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threshold is similar to the pole, D, in Fig.3. A crucial difference is that in the present paper 
the pole, which is responsible for the threshold enhancement of the spectral function in the 
sigma channel, is different from the one for the peak at zero temperature, while in ref. [7] the 
former is identihed to have continuously moved from the latter as the temperature increases. 
Which scenario is closer to reality should be clarihed in future. Also pointed out in ref. [6, 7] 
was that the range of the validity of the model is severely restricted by the unavoidable 
presence of a tachyonic pole. In our approximate solution, however, there is no tachyonic 
pole within discussed region of the parameters. Probably, the problem of a tachyonic pole 
will arise when higher loop terms are taken into account. In ref. [8], Dobado et al. studied 
the mass and the width of a and p using chiral unitary approach. They searched for the 
poles of the T-matrix of tttt scattering and showed that the width grows with temperature, 
while the mass decreases slightly. The behaviour of the pole in the a channel in ref. [8] is 
consistent with the behaviour of the pole C in Fig.3. However, they did not discuss the the 
pole, which corresponds to the pole D in Fig.3. Therefore, it is important if such a pole 
exists or not in their approach. 

B. IpI dependence 

In this subsection we discuss how the spectral functions and the poles of the propagators 
depend on the spatial momentum with the temperature hxed (T = 50, 145 and 170 MeV). 

1. TT channel 

Figs.4, 5 and 6 show the momentum dependence in the tt channel at T = 50, 145 and 
170 MeV, respectively. When T = 50MeV the pole (A) moves down at the beginning, then 
to the right and then up again as the momentum increases, approaching the real axis as 
IpI —cx). As a consequence, the peak becomes broader up to about |p| = 500MeV and then 
becomes narrower again, approaching the (5-function as |p| —> cx). At higher temperatures, 
T = 145 and 170 MeV, similar momentum dependence is observed for the spectral function 
and the pole of the propagator, but as the temperature increases the variation of the spectral 
function and the movement of the pole becomes larger. On the other hand, the pole (B) 
moves right and away from the TTr —a threshold as the momentum increase at 145 MeV. 
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As a consequence, the threshold enhancement of the spectral function decreases. At T = 170 
MeV, the pole (B), which lies on the physical sheet for |p| = 0, crosses the threshold and 
moves into the unphysical sheet as \p\ increases. Accordingly, the bound state peak becomes 
the sharp threshold enhancement and the it decreases. 

2. a channel 

Figs.7, 8 and 9 show the momentum dependence in the a channel at T = 50, 145 and 
170MeV, respectively. When the temperature is as small as 50MeV the momentum depen¬ 
dence of the spectral function is very small. At high temperature T = 145MeV, the pole (D), 
which is located near the a tttt threshold for |p| = 0, moves left and away from the thresh¬ 
old and then threshold enhancement of the spectral function decreases. At T = 170MeV, 
the pole (D), which lies on the physical sheet for \p\ = 0, crosses the threshold and moves 
into the unphysical sheet as |p| increases. Accordingly, the bound state peak becomes the 
sharp threshold enhancement and the it decreases. 

It should be noted that the numerical results presented in this section are consistent with 
the discussion in the previous section. 


V. SUMMARY 

In this paper we have studied the momentum dependence of the spectral function in the 
0(4) linear sigma model. 

We analyzed the spatial momentum dependence of the imaginary part of the self-energy. 
The temperature-dependent imaginary part of the self-energy was found to vanish in the 
large momentum limit. This is because the energy of the particles in the heat bath which 
participate in the process becomes large and therefore the Bose distribution function van¬ 
ishes. However, the temperature-dependent part of the self-energy does not completely 
vanish at large momentum; the contribution of the thermal tad pole diagrams to the real 
part of the self-energy and the temperature dependence of the mass parameter due to the 
optimization procedure do not vanish at large momentum. 

We then calculated the spectral functions and searched for the poles of the propagators. 
First, we discussed the temperature dependence for zero spatial momentum. We reproduced 
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the spectral functions in both tt and a channels in the previous work and found that the pole, 
which is responsible for the threshold enhancement of the spectral function in both tt and a 
channels, is different from the one for the peak at zero temperature. Secondly, we discussed 
the spatial-momentum dependence. When the temperature is low the spatial momentum 
dependence is also small. As the temperature becomes higher the spatial momentum depen¬ 
dence becomes also large. When the spatial momentum is smaller than the temperature, 
the spectral functions do not deviate much from those at zero spatial momentum. Once 
the spatial momentum becomes comparable with the temperature, the deviation becomes 
considerable. As the momentum becomes much larger than the temperature, thermal effects 
effectively decrease. 

These observations have an implication for heavy-ion collision experiments: When one 
extracts information about the change of intrinsic properties of hadrons at hnite temperature 
from their spectral functions measured in high energy heavy-ion collision experiments, one 
should keep in mind spatial momentum dependence discussed in the present paper. 

We expect that the essential part of the results in the present paper holds for other 
hadrons, since it is based basically on the kinematics of the particles under consideration 
and thermal particles in the heat bath. 

Finally, a comment on the higher-loop effects is in order. At finite temperature there 
exist various physical processes which are not taken into account in the present one-loop 
calculation. If we include those processes in the calculation, the structure of the spectral 
functions might be signihcantly modihed by the following reason: The sigma meson at zero 
temperature has large width due to the strong coupling with two pions. It was found in 
the one-loop calculation that the spectrum of a at hnite temperature is enhanced near the 
d —TTTT threshold. This is because as the temperature increases the mass of a decreases 
while that of tt increases. However, what is not taken into account in the one-loop calculation 
is the effect of the thermal width of tt in the process a tttt. This is taken into account in 
the calculation to two-loop order. If this effect is included this effect the sharp threshold does 
not exist any more and therefore the enhancement might be smeared out. The extension of 
the calculation to two-loop order is now in progress [12]. 
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FIG. 1: Feynman diagrams representing one-loop self-energy for tt and a. Solid and dashed lines 
correspond to a and vr respectively. 

APPENDIX: ONE-LOOP SELF-ENERGY AT FINITE TEMPERATURE 
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FIG. 5: Spectral function and the pole position of propagator at T = 145MeV for |p|=0, 145, 
1000, and \p\ ^ oo in the tt channel. 
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FIG. 8: Spectral function and the pole position of propagator at T = 145MeV for |p|=0, 145, 500 
1000, and \p\ ^ oo in the a channel. 
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